arXiv:astro-ph/0512582vl 22 Dec 2005 


Effect of Reconnection Probability on Cosmic (Super)string Network Density 

A. Avgoustidis^E and E.P.S. Shellard 1 ^ 

department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, 
University of Cambridge, Wilberforce Road, Cambridge CB3 OWA, United Kingdom 

We perform numerical simulations of cosmic string evolution with intercommuting probability P 
in the range 5 x 10~ 3 < P < 1, both in the matter and radiation eras, using a modified version of 
the Allen-Shellard code. We find that the dependence of the scaling density on P is significantly 
different than the suggested p oc P _1 form. In particular, for probabilities greater than P ~ 0.1, 
p(l/P) is approximately flat, but for P less than this value it is well-fitted by a power-law with 
exponent 0.6^o 12 • This shows that the enhancement of string densities due to a small intercom- 
muting probability is much less prominent than initially anticipated. We interpret the flat part of 
p(l/P) in terms of multiple opportunities for string reconnections during one crossing time, due to 
small-scale wiggles. We also propose a two-scale model incorporating the key physical mechanisms, 
which satisfactorily fits our results over the whole range of P covered by the simulations. 
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I. INTRODUCTION 

Cosmic strings have recently been reincarnated in the 
form of ‘cosmic superstrings’ or, more specifically, funda¬ 
mental strings and one-dinrensional D-branes produced 
at the end of brane inflation nuns (for reviews 
see for example J| 0, (gj). These objects have prop¬ 
erties which can be significantly different than usual 
field theory strings, and this opens up the possibility 
that cosmic string observations could yield information 
about string theory, while providing new ways to con¬ 
strain various brane inflation models 00 Cosmic su¬ 
perstrings typically have smaller tensions, in the range 
10- 12 < Gp < 10- 7 mum , and reconnect with proba¬ 
bilities which can be significantly less than unity [TlUl3 |. 
unlike ordinary cosmic strings 0 This reduced re¬ 
connection probability results in an enhancement of the 
string number density today maul- Cosmic super¬ 
string networks can also consist of more than one type of 
string, which can zip together to produce trilinear ver¬ 
tices with links stretching between them 110 . 

The evolution of cosmic superstring networks has re¬ 
cently been studied in Refs. 0 , 0 , 0 , 0 , 0 ]- The key 

issue is whether these networks reach a scaling regime , 
that is one in which their characteristic lengthscale stays 
constant relative to the horizon. This self-similar evo¬ 
lution is a well-known property of usual strings, which 
reconnect with probability of order unity, at least in the 
abelian case where no trilinear vertices can be formed 
(though some evidence for non-abelian scaling has also 
been found 0000). In all the above studies of 
various models of cosmic supersting networks, evidence 
for scaling behaviour was found, though the issue is not 
yet completely resolved. 

Assuming that an attractor scaling solution exists, an 
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interesting question to ask is how the density of the scal¬ 
ing network depends on the intercommuting probability 
P. Using a simple one-scale model 0, Jones, Stoica and 
Tye have argued in Ref. 0 that the expected behaviour 
is p oc P~ 2 . This is a dramatic effect, as a probabil¬ 
ity of 10 , for example, would lead to an enhancement 

of the string density by a factor of 10 4 compared to or¬ 
dinary strings: we could live in a universe packed with 
F /D-string relics from a brane inflation era, which could 
in principle be observed to extract information about 
physics at the string scale! However, string intercommut¬ 
ing is a small-scale process and can be expected to depend 
crucially on small-scale wiggles on strings, which are not 
captured by the simple one-scale model. One could still 
use a one-scale (or VOS 00) model but introduce 
an effective intercommuting probability P e g = f(P) 0 |, 
whose form should be determined by numerical simula¬ 
tions or an analytic model for small-scale structure. One 
should expect that, since small-scale wiggles would mean 
that strings may have more than one opportunity to re¬ 
connect in each encounter, the effect of P on string den¬ 
sity could be somewhat counterbalanced by the strings’ 
small-scale-structure. Thus a power law weaker than 
p oc P~ 2 may be anticipated. 


Indeed, Sakellariadou recently performed flat space 
simulations of cosmic strings with reconnection proba¬ 
bilities in the range 10 -3 < P < 1 and found a power 
law p oc P _1 in the range 10 -3 < P < 0.3 16]. The 
purpose of this letter is to investigate the dependence of 
the string energy density p on the intercommuting prob¬ 
ability P for strings evolving in more realistic cosmolog¬ 
ical backgrounds. We perform numerical simulations in 
both the matter and radiation era, and find a power law 
significantly weaker than that of Ref. 0 |. We explain 
our results in terms of string small-scale structure and 
the importance of a second scale in the problem. We 
thus propose a two-scale analytic model which provides 
a good fit to our numerical results. 
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II. SIMULATIONS 

We have performed numerical simulations of strings, 
with reduced intercommuting probabilities and evolving 
in FRW spacetime, using a modified version of the Allen- 
Shellard code 01 ■ The code sets up Vachaspati-Vilenkin 
initial conditions [27| in a horizon volume within a fixed 
comoving box, and evolves the system until the comoving 
horizon r grows to half the size of the box, when the evo¬ 
lution is stopped. The initial networks had a resolution of 
16 points per correlation length, and constant resolution 
in physical coordinates was enforced. Each network was 
given a different intercommuting probability in the range 
5 x 1CU 3 < P < 1 and was evolved both in the matter 
and radiation eras, for a dynamical range (rfi na i/ri n itiai) 
of order 3, taking several days of CPU time. 

For each network characterised by a given intercom¬ 
muting probability P, approximately ten runs were per¬ 
formed, each with different initial string density (or pre¬ 
cisely a different initial horizon to correlation length ra¬ 
tio). By plotting the time evolution of the string density 
for different initial conditions, one can bracket the scaling 
solution, as shown in Fig© getting successively more ac¬ 
curate convergence with subsequent runs. We have thus 
obtained, within errors, the scaling density p of these 
networks, each characterised by a given value of P. 

In Fig. © we plot the dimensionless parameter pt 2 /p 
(where p is the string tension) versus the inverse inter¬ 
commuting probability 1/P for matter and radiation era 
runs. In the former case P ranges from 5 x 10~ 3 to 1, but 
in the latter our limited dynamical range has at present 
only allowed us to bracket the scaling densities for P in 
the range 0.1 < P < 1. We see that for probabilities 
greater than P ~ 0.1, the function p(l/P) is approxi¬ 
mately flat, but for smaller P it develops a constant slope, 
on a log-log scale, at least in the matter era. A weighted 
fit gives a slope of 0 . 6 + 0^2 for the matter runs, and the ra¬ 
diation era data are consistent with this picture, though 
more data points are needed to confirm the value of the 
slope. Comparing to the p oc P~ 2 and p oc P _1 forms 
of Refs. 0 and JT(| respectively (also plotted), we see 
that the enhancement of string densities due to a reduced 
intercommuting probability is much less prominent than 
what was initially anticipated. For example, a probabil¬ 
ity of 5 x 10 -3 leads to an enhancement in p by only a 
factor of 10, to be contrasted with the predictions of 200 
(resp. 10 4 ) obtained from p oc P _1 (resp. p oc P -2 ). 

We have also investigated the effect of reducing the 
intercommuting probability P on the amount of string 
small-scale-structure J2S|. This can be quantified by the 
dimensionless effective tension defined as p = p e g/p, 
where p e g is the effective energy per unit length at the 
scale of the correlation length £ (this is the average length 
beyond which string directions are not correlated). For 
wiggly strings, p is greater than unity, reflecting the 
fact that string structure at sub-correlation length scales 
‘renormalises’ the tension at the scale of the correlation 
length to p e g > p (refer to Martins & EPS [0 for further 
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FIG. 1: Dimensionless string density plotted against con¬ 
formal time r for a network with intercommuting probabil¬ 
ity P = 0.75. Each curve corresponds to a different initial 
horizon to correlation length ratio. The asymptotic curves 
bracket the scaling solution, which can be estimated to be 
pt 2 /p = 3.6/°; 2 . 



FIG. 2: Dimensionless string scaling density plotted against 
the inverse intercommuting probability 1/P, for matter and 
radiation era runs. The constant slope part of the matter era 
data can be fitted by a power law with exponent OA/Iq/J- 
The overall dependence of p on P is much weaker than the 
previously suggested p oc 1/P 2 and p oc 1/P 2 forms. 

discussion of these points). In Fig. [3] we plot the dimen¬ 
sionless effective tension p against physical distance for 
a range of intercommuting probabilities. Physical scales 
are expressed in units of the physical time t. As the in¬ 
tercommuting probability decreases, we see a significant 
increase in the effective tension, especially at the scale 
of the correlation length. This is strong evidence for the 
accumulation of small-scale structure as P is reduced. 
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FIG. 3: Dimensionless effective string tension jl plotted 
against physical distance (measured in units of physical time 
t) for matter era networks with P = 1,0.1,0.05 and 0.01. 
As P is reduced, the effective tension increases signalling the 
accumulation of small-scale structure. 

III. DISCUSSION 

We can interpret the flat part of Fig. [3 in terms of 
string small-scale-structure: if two colliding strings are 
wiggly, they may have more than one opportunity to re¬ 
connect during their encounter. On large scales, correla¬ 
tion length segments move together slowly and coherently 
(v c « 0.15 for P = 1 HU), but their small-scale struc¬ 
ture is relativistic and rapidly oscillates (i> rms ~ l/\/2). 
One therefore expects the strings to cross several times 
as they pass, thus the overall process should be described 
by an effective probability of reconnection P e g > P. Re¬ 
ducing P by a factor of a few, for example, should not 
have a big impact on the string density jlj|. Indeed, we 
can understand the flat part of Fig. [3 and in particular 
explain the position where the slope starts to change, by 
using the following simple model for string interactions. 

Consider two colliding straight strings, described on 
small-scales by a monochromatic mode of amplitude A 
and period T. If v c is their relative coherent speed, the 
duration of the collision is approximately At « 2 A/v c . 
Now on small scales the strings are oscillating with period 
T, so the number of reconnection opportunities during 
their crossing time is IV ss At/T = 2 A/v c T. But since on 
these scales string modes are propagating with velocity 
c ~ 1, T is equal to the wavelength of the mode, and so 
A/T is a measure of string wiggliness. In particular, for 
large amplitudes we may write 1 + A/T ~ jl from which 
we can infer that 


Interestingly, this simple estimate gives about the right 
value for the position where the curve changes slope in 


Fig. 0 Indeed, from Fig. [3 we see that for P > 0.1, 
jl is in the range 1.6 < jl < 1.7, while we know that 
string velocities at the scale of the correlation length are 
approximately v ~ 0.15. This gives N ~ 10 intercom- 
muting opportunities. One therefore expects that to see 
a significant effect on the scaling string density, P should 
be reduced by a factor greater than about 10. 

Furthermore, we can estimate the effective probability 
by considering the N reconnection opportunities, assum¬ 
ing each has an independent probability P, that is 

Peff = 1 - (1 - P) N - (2) 

For N ~ 10, a probability as low as 0.2 yields P e s ~ 1, 
and significant deviations from unity only occur for P — 
0.1 for which P e g ~ 0.65, as observed in Fig. [3 In the 
limit P<1, Eq. (J3 gives P e g ~ NP ~ [2 (jl — l)/t; c ]P 
and so P can only be enhanced by a factor of order 10 
(see Fig. 0, which results in an effective reconnection 
probability that is still much less than unity. 

We now turn to explaining the constant slope part of 
Fig. [3 for small probabilities. First recall that a sim¬ 
ple one-scale model with P < 1 predicts p oc P~ 2 in 
direct contradiction with our numerical results. As dis¬ 
cussed above, the introduction of an effective intercom- 
muting probability can dramatically change this for rel¬ 
atively large P (say P > 0.1), thus explaining the flat 
part of Fig. [3 However, to accommodate the constant 
slope part for P < 1 one would need P e g oc P K , with 
k a constant less than 1/2. At present there seems no 
motivation for such a P e ff; instead, Eq. 0 suggests that 
P e ff is linear in P for P <1. Thus, it appears that no 
one-scale model can fit our numerical results for P 1: 
as the string intercommuting probability decreases, the 
one-scale approximation becomes increasingly poor. This 
is in agreement with the results of Ref. HU in which long 
string intercommutings were switched off in numerical 
simulations of evolving strings, while small loop produc¬ 
tion was allowed. It was found that this prevented the 
string density from scaling, that is, the inter-string dis¬ 
tance L was no longer proportional to t (the solution 
being p oc L~ 2 oc U 7 / 8 ). Nevertheless, the actual corre¬ 
lation length £ along the string did scale at approximately 
the size of the horizon (£ oc t). 

This is very similar to the situation we observe in our 
simulations. Reducing the intercommuting probability 
leads to an increase of the string density and therefore 
a decrease of the characteristic length associated with it 
(the interstring distance L). The correlation length £, 
however, is comparatively unaffected by this, and stays 
at a scale of order the horizon (Fig. 0. This can be un¬ 
derstood by considering the two distinct mechanisms for 
producing loops: (i) self-intersections of the same string 
and (ii) collisions between long strings. The former tends 
to chop off small loops, straightening the strings out (thus 
affecting mainly the correlation length £) while also con¬ 
trolling the amount of small-scale structure jl. Small loop 
production alone is a significant energ y lo ss mechanism, 
but it is not sufficient to ensure scaling [23. On the other 
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hand, long string reconnections have much more dra¬ 
matic effects introducing large-scale ‘bends’ in the strings 
which catalyse the collapse of large regions of string and 
the formation of many more loops. In contrast, these en¬ 
ergy losses are sufficient to govern the interstring distance 
L and cause scaling. 

Small loop production, the first mechanism, is not 
greatly affected by reducing the intercommuting proba¬ 
bility: once the string is sufficiently wiggly, left and right 
moving modes which fail to interact due to a small P 
will keep propagating and meet more incoming modes, 
with which they will eventually interact and form loops. 
It might take longer, but such interactions are inevitable 
even for P< 1. The relevant question is whether enough 
of these interactions can take place in each Hubble time 
in order to straighten the strings out at the horizon 
scale, but this seems reasonable given the much shorter 
timescale on which small loop production operates. It 
appears to be confirmed by Fig. 0 where there is only 
a very weak build-up of small scale structure (a factor 
of 2-3 in p) while P changes by over two orders of mag¬ 
nitude. On the other hand, long string intercommuting, 
the second mechanism, depends more crucially on P: two 
colliding strings have a given interval of time in which to 
interact, so intercommuting is no longer inevitable for 
relatively small P < 0.1. A reduced P necessarily means 
less string interactions and less energy dumped from the 
long string network in the form of loops. This leads to 
an increase in the long string density and thus a signifi¬ 
cant decrease of the characteristic length relative to the 
correlation/horizon scale. 

To obtain an analytic model for such networks, it is 
therefore necessary to introduce two scales: a character¬ 
istic length L (roughly the interstring distance) quantify¬ 
ing the energy density in strings, and a correlation length 
£, defined to be the distance beyond which string direc¬ 
tions are not correlated. Indeed, similar models have 
appeared in the literature (e.g. the three-scale model of 
Ref. 0), but for normal P = 1 cosmic strings it has been 
established that these two scales are comparable so the 
one-scale (velocity-dependent) approximation mum 
can be surprisingly accurate [32]. Here, we develop a two- 
scale velocity-dependent model which fits our numerical 
results, again surprisingly well given its simplicity. 

Consider a string network, characterised by tension p, 
correlation length £ and energy density p. We define the 
characteristic lengthscale L of the network by p = p/L 2 
and note that the number of strings per correlation vol¬ 
ume V = £ 3 is = pV/pt; = £ 2 /P 2 . If v is the typical 
string velocity, each string intersects N^ — l other strings 
in time At = £/u, so we have _/V| — Ng intercommut- 
ings per correlation volume per time At. Assuming that 
each intercommuting produces a loop of length c£ (this 
can be done formally by integrating an appropriate loop 
production function over all relevant loop sizes, which 
introduces the loop production parameter c [HI) the en¬ 
ergy loss due to the formation of loops can be written 


as 


(-) 

\ 5t J 2—scale 


(N 2 - N 6 )v 
£ 4 


/xc£ = cp 




(3) 


This is to be contrasted with the corresponding result for 
the one-scale model 


(*£) 

\^t J x_ sca i e 



(4) 


As an interesting aside, we note that the last term in 
Eq. 0 oc l/£ has the same form (though opposite sign) 
as the term which would need to be introduced to account 
for direct small loop production (or string radiation); it 
cannot itself cause L to scale as t. 

Our velocity dependent two-scale model can therefore 
be constructed by using the usual VOS model equations 
[25j , derived by performing a statistical averaging proce¬ 
dure on the Nambu-Goto equations of motion and energy 
momentum tensor, but using the phenomenological loop 
production term m instead of ©• The result is a sys¬ 
tem of two coupled ODEs, governing the time evolution 
of the characteristic length L and the average velocity v 
of string segments: 

„g?P _ /£ L\ , . 

= 2 a Hl + V ) + CV \L~y (5) 

dv , n. f k h \ 

1 , = (1 - *’■> u - Vj (6) 

In Eq. k is the so-called momentum parameter, which 
is a measure of the angle between the curvature vector 
and the velocity of string segments and thus is related 
to the smoothness of strings [24j. Note that we have 
made no attempt to derive an evolution equation for the 
correlation length £ as the simulations show that this 
depends weakly on P and moreover, unlike L , it remains 
comparable to the horizon (Fig. 0. We will neglect this 
weak dependence of £ on P and take £ = t in Eqs. 00 . 

To apply the model 00 to our simulated networks, we 
have also introduced an effective intercommuting proba¬ 
bility given by ©, that is, we have replaced c in Eq. © 
by P e fiC. As a first approximation, we have neglected the 
dependence of N on P (this can be understood in terms 
of the dependence of p on P, see Eq. 0 and Fig. 0 and 
we have simply taken N = 10 (a rough average value over 
the range of P). For the loop production parameter, we 
have used the value c = 0.23 of Refs. fpj. I.32t which fits 
both radiation and matter era runs in the P = 1 case. 

In Fig. 0 we plot the matter era scaling energy den¬ 
sities, obtained both from the simulations and our two- 
scale model (solid line), for string networks with inter¬ 
commuting probabilities in the range 5 x 10” 3 < P < 1. 
We see that the model provides a surprisingly good fit 
(given its simplicity and the approximations made) of the 
numerical data and it reproduces the observed change of 
slope around P ft 0.1. We stress that we have made 
no special parameter choices to obtain this fit, and by 
modifying parameters further we could do much better. 
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Nevertheless, it is clear from the dashed line in Fig. Q] 
how the model could be improved by taking into account 
the dependence of the collision number N on p{P). In 
this case, we have estimated the effective string tension 
p as a function of P from Fig. [3 which yields N from 
Eq. (JTJ (that is, 6-18 over the full range of P) and then 
P e g(P) from Eq. © which is finally input ‘by hand’ in 
our two-scale model. The striking agreement that re¬ 
sults provides strong motivation for understanding the 
key physical mechanisms governing the small-scale struc¬ 
ture parameter p. Indeed, modifying P seems to provide 
a useful testbed for observing the dynamics of p, but we 
leave a more sophisticated analysis for a future publica¬ 
tion §14} • We note the need also to take into account the 
dependence of £ on the reconnection probability, and a 
direct small-loop production term (as discussed above). 
However, here, we just highlight the fact that a simple 
version of our two-scale model with the addition of the 
effective probability of Eq. 0 seems to provide a satis¬ 
factory fit to the numerical data. 



FIG. 4: Scaling string density obtained from simulations 
(data points with errors) and from the analytic two-scale 
model irilfil (solid line). The fit is improved further (dashed 
line) by phenomenologically incorporating the dependence of 
the effective string tension p on the reconnection probability 
P. 


scale) of about 0.6lg^2- This yields very different pre¬ 
dictions from the generally expected p oc P form; a 
probability of P = 5 x 10 -3 , for example, gives an en¬ 
ergy density enhancement of a factor of 10, compared 
to a factor of 200 with p oc P _1 . Clearly, the results 
presented here are important for determining the quan¬ 
titative observational predictions of models with cosmic 
(super-)strings. The distinction between ordinary cosmic 
strings and strings with P < 1 is obviously more subtle 
than first anticipated. In particular, we note that cosmic 
(super-)strings will not be as strongly constrained ob- 
servationally and many limits in the literature will have 
to be re-examined in light of these results (or, at least, 
re-normalised). 

We have also endeavoured to provide some physical 
explanations for why the string density depends non- 
trivially on the reconnection probability. We can explain 
the flat dependence of p(l/P) in terms of small-scale wig¬ 
gles on strings, which lead to multiple opportunities for 
reconnection in each crossing time, thus introducing an 
effective reconnection probability P e g = f(P). By ap¬ 
proximating the small-scale structure on long strings by 
a monochromatic mode, we suggested a physically moti¬ 
vated form for P e ff, which adequately explains the large 
P behaviour of p(l/P). 

The simulations demonstrated that, although the 
string density increases for small P leading to a reduction 
of the characteristic length scale (or interstring distance 
L) of the network, the correlation length £ (the length be¬ 
yond which string directions are not correlated) scales at 
a size comparable to the horizon. We have explained how 
one can understand this fact in terms of the two distinct 
mechanisms for loop production: the self-intersection of 
wiggly strings which produces small loops and tends to 
straighten the strings, determining the correlation length 
£, and long-string collisions, which cause much greater 
loop energy losses, determining the characteristic length 
L. We used a simple two-scale model to describe the 
behaviour of p(l/P) over the whole range of probabili¬ 
ties probed by the simulations. Our model provides an 
adequate fit to the numerical data, given its simplicity 
and the approximations used. We leave a more detailed 
numerical investigation and further improvements of the 
analytic model for a future publication. 


IV. CONCLUSION 

We have performed numerical simulations of cosmic 
strings with intercommuting probabilities in the range 
5 x HT 3 < P < 1 evolving in a matter- or radiation- 
dominated FRW universe. We have found that the de¬ 
pendence of the string density p on the intercommuting 
probability P is much weaker than previously suggested 
in the literature. In particular, the function p(l/P) has 
an initially flat dependence, up to probabilities of about 
0.1, and then develops a constant slope (on a log-log 
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